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calculation is implemented in the Monte Carlo event generator HORACE, which is available 
for data analysis. 
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1. Introduction 

The Drell-Yan-like production of single W and Z bosons, with the weak boson decaying 
into a lepton pair, is a clean process with a large cross section at hadron colliders. It is 
well suited for a number of precision measurements both at the proton-antiproton (pp) 
Fermilab Tevatron Run II and at the proton-proton (pp) CERN Large Hadron Collider 
(LHC) While the charged-current process pp(pp) — > W — > lv\ + X (I = e, /i) 

can be used for a precise determination of the VF-boson mass and width, the neutral- 
current process pp(pp) — ► 7, Z — > / + /~ +X is of interest because Z data help to accurately 
calibrate detector parameters, such as the energy scale and resolution of the electromagnetic 
calorimeter, which are essential ingredients for a precise determination of the W mass 
^, |3|, The neutral-current process can also be used to determine the W mass from 
the ratio of the transverse mass distributions of the W and Z boson |J (especially at 
high luminosity), as well as to extract the effective weak mixing angle from the forward- 
backward asymmetry ||. Furthermore, both processes are backgrounds to many new 
physics searches and can also be useful to monitor the collider luminosity and measure the 
Parton Distribution Functions @. 

The success of precision studies in W and Z production at the Tevatron and the LHC 
depends on the progress in reducing theoretical uncertainties, if possible to the 1% level. 
To this end, QCD and electroweak corrections must be under control. Recent developments 
in the theory of W and Z production at hadron colliders are reviewed in Ref. ||. Order 
a QED corrections to the leptonic decays of W and Z bosons were calculated in Ref. B. 
For the neutral-current process, which is the concern of the present study, the complete 



calculation of 0(a) QED corrections was presented in Ref. [lCj, while the full set of one-loop 
electroweak corrections was computed in Ref. [p]] . To evaluate the impact of electroweak 
corrections in the experimental analysis, the CDF and D0 collaborations at the Tevatron 
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Run I made use of the fixed-order calculations of Refs. |], [T^] for single W production 1 
and of Refs. g, 0, 0] for the single Z process. However, the anticipated precision for 
the Drell-Yan process at the Tevatron Run II and the LHC and, in particular, for the 
high-precision measurement of the W mass, which depends on a precise determination of 
the Z parameters, requires that leading contributions from radiation of multiple soft and 
collinear photons are resummed through all orders of a |8], 10 1. A first attempt towards 
the inclusion of higher-order QED corrections was the calculation of W and Z production 
with radiation of two additional photons in Ref. 14]. Recently, higher-order corrections 
due to multi-photon (real and virtual) radiation in W decays at hadron colliders have 
been computed independently in Refs. fl5| , 16|. Previous investigations and Monte Carlo 
simulations of multiple QED radiation in W/Z production in hadron collisions can be found 
in Ref. [17]. A first step towards the inclusion of QED and QCD radiation into a single 
generator has been recently performed in Ref. by combining QCD resummation with 
0(a) final-state QED corrections. 

The aim of the present work is to compute higher-order final-state QED corrections 
to the neutral-current Drell-Yan process pp(pp) — > 7, Z — > l + l~ + X and, in particular, 
to quantify the effects induced by such corrections on the experimental observables of in- 
terest in view of a precise determination of the W mass at hadron colliders. We restrict 
ourselves to consider final-state QED corrections only, because it is known from previous 
investigations that electroweak corrections to Drell-Yan-like processes are largely domi- 
nated by photon radiation from the final-state charged leptons [1C, 11, [l^, [l3|], because of 
the presence of large lepton-mass logarithms arising from collinear radiation. 

The present hadron collider average for the mass of the W boson M\y is Myy = 
80.454 ± 0.059 GeV £j]|. The target M w precision is of the order of 30 MeV (per channel 
per experiment) for the Tevatron Run II and of about 15 MeV for the LHC |l|, ^], which 
corresponds to a high-precision determination of M\y with a relative accuracy at the 10~ 4 
level. The shift in the Z mass due to multi-photon radiation was estimated at the Tevatron 
Run I to be less than 10 MeV by the CDF collaboration for the Z — > decay Q (when 

using the package PHOTOS of Ref. @) and 10 MeV by the DO collaboration for the 
Z — > e + e~ decay f|]. These shifts are presently treated as systematic uncertainties on the 
Z mass. However, in view of the expected precision for M\\r in Run II and at the LHC, it 
will be necessary to include multiple photon corrections into account when extracting the 
W mass from data or when calibrating detector components using Z data. The present 
paper aims at completing the recent calculations of higher-order QED corrections to single 
W production [jl5|, |l6| to cover the process of single Z production, in order to eliminate 
(or largely reduce) this source of theoretical uncertainty in the experimental analysis. 

The paper is organized as follows. In Sect. ^ we describe our theoretical approach. In 
Sect. |3| we present the phenomenological implications of our study, by discussing compar- 
isons with available calculations and quantifying the effects of (higher-order) QED correc- 
tions to a number of distributions and observables of experimental interest. In Sect. |I] we 



1 The calculation of Ref. [jl2j refers to electroweak corrections contributing in the so-called pole approx- 
imation. Complete calculations of O(a) electroweak corrections, beyond the pole approximation, to single 
W production are given in Ref. (13]. 
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give our conclusions and perspectives. 



2. Theoretical approach 

The observable cross section for the Drell-Yan process at haelron colliders with nominal 
squared centre of mass (cm.) energy s is obtained by convoluting the cross section of 
the parton process qq — > 7, Z — > I = e,fi, with the Parton Distribution Functions 

(PDFs) and summing over all quark flavors q. According to factorization theorems, it can 
be written as 2 

°"( s ) = / dcos9dxidx 2 dx 3 dx 4 (f q / A (xi,Q 2 ) fq/ B (x 2 ,Q 2 ) + (q <-> q)) 

q J 

Dl -(x 3 ,Q 2 )D l+ (x 4 ,Q 2 ) d ^° s9 f (2.1) 

a cos 

where da (cos 6, s)/d cos 9 is the parton-level differential cross section, function of the scat- 
tering angle 9 in the cm. frame and the squared parton cm. energy s = xix 2 s. The 
function f q {q)/{A,B)( x ii Q 2 ) stands for the PDFs of the initial-state quarks (antiquarks) with 
momentum fractions Xi (i = 1, 2; < xi < 1) inside the proton(antiproton), where (A, B) = 
(p,p) for the Tevatron and (A,B) = (p,p) for the LHC. The quantities Di-^+^(xi,Q 2 ) rep- 
resent the QED Structure Functions (SF) for the final-state lepton(antilepton), which can 
be interpreted as the probability density of finding inside a parent lepton a lepton with 
momentum fraction Xi (i = 3,4) at a virtuality scale Q 2 after photon radiation The 
integration is over the parton momentum fractions x\ and x 2 and the lepton momentum 
fractions X3 and x±. This amounts to neglect the contribution of photon radiation off the 
initial-state quarks, as well as the contribution of initial-final-state photon interference. 
The inclusion in eq. ( p.l| ) of final-state photon radiation only is meaningful from the point 
of view of quantum field theory because final-state QED corrections to a neutral-current 
2 — > 2 process form a gauge-invariant subset within the full set of electroweak corrections 
and can be therefore treated separately. Furthermore, photon radiation off (initial-state) 
quarks gives rise to quark mass singularities which must be reabsorbed in PDFs, in analogy 
to gluon emission in QCD [p3| 3 . After this "renormalization" procedure, it is known from 



previous investigations |10| , 11, 12, 13] that 0(a) QED initial-state and initial- final-state 
interference corrections to the relevant distributions are uniform and very small (at the 
0.1% level), while final-state photon corrections are seen to completely dominate and are 
responsible for strong modifications, at the level of several %. Consequently, initial-state 
and initial-final-state interference corrections contribute very little, for instance, to Z-boson 
mass shifts, as demonstrated in Ref. [pX|] , and can be safely neglected for the purposes of 
the present study. 



2 According to eq. ([2.l|), the transverse motion of the Z boson is neglected. The modeling of Z transverse 
momentum requires careful QCD calculations, including the resummation of multiple soft-gluon radiation, 
as available in the Monte Carlo program RESBOS Bl used at the Tevatron. 

3 Very recently, the MRST group performed a global analysis of PDFs including QED corrections and con- 
sistently incorporated in the parameterization of PDFs the contribution of QED radiation from quarks [pij . 
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To compute higher-order photon corrections and simulate the mechanism of multi- 
photon emission, we make use of the QED Parton Shower (PS) approach plfl . It consists 
in a numerical solution of the QED Gribov-Lipatov-Altarelli-Parisi evolution equation for 
the lepton SF D(x, Q 2 ) in the non-singlet channel. The solution can be cast in the form pa 

D(x, Q 2 ) = U(Q 2 , m 2 )5(l - x) 

/ n \ f® 2 Ht' f x + 

+ (^J J a U(Q 2 ,s')—U(s',m 2 ) J dyP(y)S(x-y) 

+ | — ] / 11(C)-.,-')— I H(,'., ")_H(,", m 2 ) x 



z,n / Jm 2 b Jm 2 b 



■ + 



dx! / dx 2 P{ Xl )P(x 2 )5{x - Xl x 2 ) + • • • (2.2) 



where P(x) 



1 4- r 2 f l \ + z 2 

P(x) = - 5(1 -x) dz (2.3) 

1- x J 1 - z 



is the regularized splitting function for the QED branching I — > Z7 and the factor II 

a f Sl ds' f x + 
2n 



n(si,s 2 ) = exp 



dzP(z) 



(2.4) 



s 2 b JO 

is the Sudakov form factor, representing the probability that a lepton evolves from virtu- 
ally s 2 to virtuality s% with no emission of photons of energy fraction greater than (an 
infrared regulator) e = 1 — x + . Equation (2.2) allows to compute D(x,Q 2 ) by means of a 
Monte Carlo algorithm which, as shown in detail in Ref. p5| , simulates the emission of a 
shower of (real and virtual) photons by a lepton and accounts for exponentiation of soft 
photons and resummation of collinear logarithms due to multiple hard bremsstrahlung. An 
0(a) expansion of the complete PS algorithm can be worked out as well, as described in 
detail in Ref. |25|], in order to compare with fixed-order calculations. A clear advantage of 
the PS algorithm with respect to a strictly collinear approximation is the possibility of gen- 
erating transverse momentum px of leptons and photons at each branching. In the present 
implementation of the PS algorithm, the generation of the photons' angles is performed 
according to the factorized part of the Z radiative decay matrix element Z — > / + / _ 7, i.e. 

2 



cos 7OC _(_£^_ (2.5) 



Pl+ Pi_ 

where Pi-n+\ is the lepton(antilepton) four-momentum and k is the photon four-momentum. 
The generation of transverse degrees of freedom at each branching allows an exclusive event 
generation suitable to implement experimental cuts according to a realistic event selection, 
as shown in the next section. 



3. Numerical results and discussion 

The formulation described above has been implemented in the MC event generator HO- 
RACE, developed in Ref. [16| to quantify the effect of higher-order final-state QED correc- 
tions on the VF-mass determination at hadron colliders. The predictions of HORACE were 
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Table 1: The fraction of events (in %) with a photon of energy greater than E™ n = fco-Ebeam 
as predicted by the present MC program (HORACE) and the calculation of Ref. j| (B&K), at a 
parton cm. energy of 90 GeV. 



recently compared with those of the independent MC generator WINHAC [15] for various 
single- W observables of experimental interest, finding a good agreement ||26|| . However, in 
order to check the implementation in HORACE of QED corrections to the neutral-current 
process, we performed further numerical tests in comparison with available fixed-order 
calculations. 



3.1 Comparisons with fixed-order calculations 



First, we compared the results for the 0(a) hard-bremsstrahlung correction to the parton 
level cross section of Ref. || with the corresponding 0(a) predictions of HORACE. The 
results of this comparison are given in Tab. [I], which shows the fraction of events (in %) 
as a function of the lower cut-off k (k = E^/E\j eam , where E y is the photon energy and 
-^beam is the parton beam energy), for both electrons and muons at a parton cm. energy 
of 90 GeV. In Tab. |] the results of Berends and Kleiss calculation of Ref. H are denoted 
as B&K. The (absolute) differences between the predictions of the two calculations for the 
fraction of hard-photon events are at a few per mille level, as expected on the basis of the 
leading logarithmic approximation inherent the PS algorithm. 

We also performed tuned tests between the predictions of HORACE and those of the 



program ZGRAD/ZGRAD2 (denoted in the following as ZGRAD2) |,§, in order to 



compare the results for hadron-level processes. The input parameters used in the compar- 
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Tevatron Run II 


LHC 


Program 


e /i 


e fjL 


ZGRAD2 EBA 
HORACE EBA 
ZGRAD2 
HORACE 0(a) 
HORACE exponentiated 


44.11(1) 44.11(1) 
44.07(1) 44.07(1) 
40.84(4) 42.42(4) 
40.65(3) 42.22(1) 
40.60(1) 42.21(1) 


557.6(1) 557.6(1) 
557.2(1) 557.2(1) 
510(1) 533(1) 
507.2(1) 530.4(1) 
507.6(1) 530.5(1) 



Table 2: Comparison between the present calculation (HORACE) and ZGRAD /ZGRAD2 |Io[ [llj 
(denoted as ZGRAD2) for the pp,pp — > 7, Z — > I = e,/j, cross sections (in pb), at the cm. 

energies of the Tevatron Run II (y/s — 2 TeV) and the LHC (y/s = 14 TeV), according to the input 
parameters and cuts discussed in the text. 



isons with ZGRAD2, as well as in the simulations described below, are the following: 

m e = 0.511 x 10~ 3 GeV m M = 0.10565836 GeV 

a' 1 = 137.03599976 G„ = 1.16639 x 10" 5 GeV -2 

M z = 91.1876 GeV sin 2 6$ = 0.23150 T z = 2.4952 GeV (3.1) 



according to the latest PDG values [27]. In terms of these parameters, we implement in the 
amplitude for the parton processes qq — > 7, Z — > l + l~ , with q = u,d, s, c, b, the running of 
the electromagnetic coupling constant, as due to fermion loop contributions to the photon 
vacuum polarization 4 , the running Z width in the Z propagator and effective vector 
and axial- vector coupling constants, in agreement with the Effective Born Approximation 
(EBA) described in Ref. |2], [ll]], which yields a good description of the complete 0{a) 
electroweak calculation in the region around the Z pole. 

The numerical results are obtained by imposing the following transverse momentum 
(j>t) and pseudo-rapidity {rf) cuts on the final-state leptons: 

PtQ?) > 25 GeV |r/(Z ± )| < 1.2 (for the Tevatron) 

Pt(^) > 25 GeV |r ? (Z ± )| < 2.4 (for the LHC) (3.2) 

which model the acceptance cuts used by the experimental collaborations in the analysis 
of the Drell-Yan process. In the comparison with ZGRAD2, we impose an additional cut 
on the lepton-pair invariant mass M ; +;-, i.e. 75 GeV < < 105 GeV, to isolate 

the energy region around the Z resonance. The results refer to final-state leptons in the 
absence of lepton identification requirements, i.e. for so-called bare leptons. 

The predictions for the hadron-level processes are obtained by convoluting the parton- 



level matrix element with the CTEQ6 PDF set [29]. The scale Q 2 is set to be Q 2 = s, s in 
both PDFs and lepton SFs. The cm. energies considered are y/s = 2 TeV for the Tevatron 
Run II and xfs = U TeV for the LHC. 



The non-perturbative hadronic contributions to photon vacuum polarization is included in terms of the 

(5) 

hadrons 



updated parameterization of Aa[ 5 J drons of Ref. 
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Figure 1: The invariant mass distribution M[+[- (left panel) and relative effect of 0(a) (right panel, 
up) and higher-order QED corrections (right panel, down) for the Z lepton decays Z — ► e + e~, 
at y/s — 2 TeV, according to the lepton identification criteria discussed in the text. 



The results of the hadron-level comparison between HORACE and ZGRAD2 are shown 
in Tab. ^, where the first two lines are the predictions of the programs corresponding to the 
cross section in the EBA, while the third and fourth lines refer to the QED corrected cross 
section in the presence of O(a) final-state photon corrections. For completeness, we show 
in the last line also the results of HORACE when including multiple QED radiation. From 
Tab. §, it can be seen that first-order QED corrections lower the cross section of about 
10% for electrons and of about 5% for muons, both at the Tevatron and the LHC. The 
effect of multiple photon radiation on the integrated cross section is at the 0.1% level. The 
predictions of the two programs for the EBA cross sections are in very good agreement, 
being the relative differences below the 0.1%. Concerning the 0{a) QED corrected cross 
sections, the results of the two programs differ at the 0.5% level. As in the case of the 
parton-level comparison shown in Tab. this discrepancy can be ascribed to next-to- 
leading-order contributions which are included in the full 0{a) calculation implemented in 
ZGRAD2 and are missing in the PS algorithm of HORACE. 

3.2 Distributions 

Having established the physical and technical accuracy of the implementation of the neutral- 
current process in HORACE, we move to the presentation and discussion of the impact of 
higher-order corrections, with particular emphasis on the effects to the observables of in- 
terest for the measurement of the W mass, i.e. the lepton-pair invariant mass distribution 
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Figure 2: The Z transverse mass distribution M§ (left panel) and relative effect of 0(a) (right 
panel, up) and higher-order QED corrections (right panel, down) for the Z lepton decays Z — > 
e + e~,/i + /i~ at t/s — 2 TeV, according to the lepton identification criteria discussed in the text. 



Mi+i~ and the Z transverse mass distribution M| g § g |] , which are defined as follows 



V(^++^-) 2 -(p i+ +Pi-) 2 



2p T (Z+)p T (Z-)(l-cos 



(3.3) 



In eq. ( |3.3[ ) Ei~n+\ and P/-(;+) are the lepton(antilepton) energy and three- momentum, 
Pt(1~),Pt(1 + ) are the transverse momenta of the lepton and antilepton and 



d+l- 



is the 



angle between the two leptons in the transverse plane. In order to perform a more realistic 
phenomenological analysis and study the dependence of the QED corrections from detec- 
tor effects, we implement, in addition to the cuts of eq. (|3.2|), the lepton identification 



requirements quoted in Table I of Ref. [12|. According to these criteria, electron and pho- 
ton four-momenta are combined for small opening angles between the two particles, while 
muons are identified as hits in the muon chambers with an associated track consistent 
with a minimum ionizing particle. The four-momenta recombination for the electrons cor- 
responds to a calorimetric particle identification, while the identification requirements for 
muons are close to a bare event selection. All the numerical results shown in the following 
refer to the cuts of eq. ( |3.2| ) for the Tevatron, at y/~s = 2 TeV. 

The results of our MC simulations are shown in Fig. |] (invariant mass) and Fig. || 
(transverse mass), when considering for the final-state electrons and muons the lepton 
identification criteria discussed above. For the sake of comparison, we also included the 
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order a 




Figure 3: The Z rapidity distribution (left panel) and effect of O(a) (right panel, up) and higher- 



order QED corrections (right panel, down) for the Z lepton decays Z 
according to the lepton identification criteria discussed in the text. 



e+e~ , fi+n~ at y/s = 2 TeV, 



results corresponding to bare electrons. The ratios shown in Fig. |l| and in Fig. Q are 
defined as <7j, a / V^eba anci c^h.o./^cu where ct^ebAj <?i,a anci 0"i,h.o. are the MC predictions 
for the differential cross section at the EBA, 0(a) and higher-order level, for bin i, respec- 
tively. First of all, we notice that our results for the shape and size of 0(a) corrections 
to the invariant mass distribution are in agreement with the results of the diagrammatic 
calculation of Ref. ]lC|| . For both the invariant mass and transverse mass distribution 
and for bare leptons, the higher order corrections change sign when passing from below to 
above the Z peak, and the relative correction reaches the order of 10% for bare electrons 
and of 1% for muons, because the collinear logarithm am(s/m?) is larger in the electron 
case. After photon recombination, the higher-order correction for calorimetric electrons 
becomes flat and is reduced well below the 1% level, because of the disappearance of the 
lepton-mass logarithms. We also investigated the impact of higher-order corrections on 
the invariant mass distribution for large values of M1+1-, i.e. Mi+i~ > 200 GeV, which 
is important in the search for new-physics and where first-order electroweak corrections 
amount to several per cent [11|. We found that higher-order QED corrections are at a few 
% level and therefore are small in comparison with other theoretical contributions, such as 
Sudakov-like electroweak logarithms [11], and the expected statistical uncertainty in this 
invariant mass range. 

In our study, we considered other observables of interest for precision electroweak 
measurements, such as the forward-backward asymmetry [||, ||, 10, O], and the quantities 
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Figure 4: The forward-backward asymmetry Afb (left panel) and effect of 0(a) (right panel, up) 
and higher-order QED corrections (right panel, down) for the Z lepton decays Z — > e + e~,/.t + /i~ at 
^/s = 2 TeV, according to the lepton identification criteria discussed in the text. 



relevant for luminosity monitoring, e.g. the lepton and Z rapidity distributions The 
results for the Z rapidity are shown in Fig. |3| We observe that 0(a) QED corrections 
are at the level of 3-4% for muons and of 1-2% for calorimetric electrons. They are of the 
same order of next-next-to-leading order (NNLO) QCD contributions recently computed 
in Ref. [31]. The contribution of higher-order corrections is quite small, at the 0.1% level. 
This can be understood as follows. The luminosity observables are not strongly varying 
distributions and this implies that 0{a) QED corrections, and consequently higher-order 
corrections too, are almost flat and smaller with respect to the case of the invariant mass 
and transverse mass distributions. These results are in agreement with what observed 
for the analogous single- W observables in Ref. [26]. The results for the forward-backward 
asymmetry Afb at the Tevatron Run II are shown in Fig. |] as a function of the lepton-pair 
invariant mass, when using for Afb the definition for pp collisions given in Ref. p0|] . The 
effects of QED corrections shown in the right panel of Fig. B are defined as A FB - A F B B A 
and A^p B ' — A FB , giving the (absolute) contribution of 0{a) and higher-order corrections, 
respectively. Our results for the 0(a) corrections well agree with those presented in Ref. 
|1C], confirming that 0(a) QED corrections are large in the region below the Z peak and 
small around and above it. Higher-order corrections are about a factor of ten smaller, 
modifying the asymmetry of about 0.01 below the peak and 0.001 above it. 
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Table 3: The Z mass shifts due to 0(a) (AMf) and higher-order QED corrections (AM| ° ), at 
\fs = 2 TeV, according to the different experimental conditions discussed in the text. The numbers 
in parentheses are statistical errors for the last digits. 



3.3 Z-boson mass shifts 

To complete the phenomenological analysis, we investigated the shift induced by multiple 
QED radiation on the Z fitted mass, which is an important parameter for calibration 
purposes. The strategy followed by the CDF and D0 collaborations to extract Mz from 
the data is to perform a fit to the invariant mass distribution or to the Z transverse 
mass distribution, as a consistency check || In our study, we perform binned \ 2 fits 
to the Mi+i~ distribution, in analogy with the experimental procedure and following the 
approach already used in Ref. JO]] for the W mass. According to this procedure, we 
generate a sample of "pseudo-data" for the invariant mass distribution, calculating the 
distribution at the EBA level in terms of an input Z mass eba, input ^ w ith high numerical 
precision. Next, we compute the Mj+j- distribution including 0(a) radiative corrections 
for a number of Z mass values and we calculate, for each Mz value, the \ 2 as 

X 2 = " ^>eba) 2 /(A^ q + Aa? EBA ) (3.4) 

i 

where a; eba and (J% Qi 3X6 the Monte Carlo predictions at the EBA and 0(a) level, re- 
spectively, for bin i and Act^ebA; the corresponding statistical errors due to nu- 
merical integration. From the minimum of the x 2 distribution, we derive the fitted Z 
mass value M^' fitted and we quantify the mass shift due to 0(a) corrections as AM^ = 
^eba, input _ j^j-a .fitted ^ rj,^ g^fj- j ue ^ higher-order corrections is derived according 

to the same procedure, by generating a sample of "pseudo-data" for the Mj+j- distri- 
bution at 0(a) with an input Z mass value M^' mput and fitting them in terms of tem- 
plates of the invariant mass distribution obtained by including higher-order corrections for 
different Mz values. The mass shift induced by higher-order corrections is obtained as 
AM|'°' = M°;' mput — My '' tted , where ^jj-°'> fitted j s ^ ne % mass value returned by the fit. 

The results of these MC experiments are given in Tab. [3|, showing the Z mass shifts due 
to 0(a) (AM|) and higher-order (AM|°') corrections for electrons and muons detected 
according to different combinations of experimental conditions at y/s = 2 TeV. Indeed, 
in addition to the shifts obtained in the presence of the lepton identification requirements 
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discussed above and denoted as Lepton ID in Tab. |j| we also quote the results obtained 
when taking into account the uncertainties in the energy and momentum measurements 
of the leptons in the detector. These uncertainties are simulated by means of a Gaussian 
smearing of the particle four-momenta, using as standard deviation values the specifications 
relative to electrons and muons for the Run II D0 detector They are denoted as 

Smearing in Tab. ||. The different combinations of smearing effects and lepton identification 
criteria shown in Tab. |3]aim at emphasizing the dependence of the Z mass shift on detector 
specifications. In particular, it can be seen that the smearing effects tend to enhance the 
shifts, whereas the lepton identification requirements reduce them. When considering a 
specific experimental condition and comparing the 0(a) mass shift with the higher-order 
one, it can be noticed that higher-order corrections slightly reduce the effect due to O(o) 
contributions, as expected, and that the absolute value of the Z mass shift due to multiple 
photon radiation is, in the most realistic case, of the order of 10% of that caused by one 
photon emission 5 . This agrees with the conclusions of Ref. [16] for the W mass shifts due 
to QED effects. When using the same cuts and (simplified) detector model, we checked that 
the above conclusions are valid for the LHC too, being the relative effect of QED corrections 
on the invariant mass distribution of the same order as the shifts at the Tevatron. A more 
realistic analysis would require a full detector simulation, as well as a modeling of the 
transverse motion of the Z, which are beyond the scope of the present paper. 



4. Conclusions and perspectives 

Drell-Yan production of Z bosons is an important process for precision measurements at 
hadron colliders. In particular, a precise extraction of the W mass in a hadron collider 
environment requires a simultaneous determination of the Z mass in lepton pair production 
for calibration purposes. In order to perform such precision studies at the Tevatron Run II 
and the LHC, radiative corrections, including the contribution of higher-order effects, must 
be under control. In this paper, we have presented a calculation of multiple photon final- 
state corrections to the neutral-current Drell-Yan process, with the aim of contributing to 
the progress in reducing present theoretical uncertainties. The calculation is implemented 
in the MC program HORACE, which can be used for data analysis. The program has been 
carefully tested against independent fixed-order calculations, finding good agreement. 

We discussed the phenomenological implications of our calculation, concentrating on 
the observables which are important for a precise determination of the electroweak param- 
eters and for luminosity monitoring. In the presence of lepton identification requirements, 
multiple photon corrections to the lepton-pair invariant mass and to the Z transverse mass 
distribution were found to modify the distributions around the Z peak by about 0.1-1%. 
For the forward-backward asymmetry, higher-order QED corrections modify this observ- 
able by about 0.01 below the peak and 0.001 around and above it. The quantities of 

5 The opposite sign between O(a) and higher-order mass shifts can be understood by noticing that, in 
the vicinity of the Z pole, multiple photon corrections to the invariant mass distribution damp the negative 
effect of 0(a) corrections (see Fig. introducing a positive contribution which makes the exponentiated 
QED corrected distribution closer to the EBA distribution than the 0(a) one. 
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interest for luminosity, such as the lepton and Z rapidity distributions, are slightly affected 
by multiple QED radiation, at the 0.1% level. We also investigated the shift on the Z 
mass, an important parameter to calibrate detector components. The precise value of the 
shift was found to significantly depend on detector effects (as expected), being the abso- 
lute value of the mass shift due to multiple radiation of the order of 10% of that induced 
by 0(a) corrections. When associated to the conclusions of Refs. [15, [L6|, for single 
W production, the results of the present paper indicate that the effects of multiple QED 
radiation are non-negligible in view of the expected precision at the Tevatron Run II and 
the LHC and should be carefully considered in future experimental analyses. 

Concerning possible developments, it would be interesting to compare the predictions 
of HORACE with those of other programs at the level of differential distributions and, 
in particular, with the photon distributions of the package PHOTOS, which is presently 
used at the Tevatron to estimate the impact of two-photon radiation. A second possible 
perspective would be the merging of exact 0(a) electroweak corrections with the exclusive 
photon exponentiation of the PS algorithm, to improve the physical precision of present 
calculations. A further, long-term development is the combination of electroweak and QCD 
radiation into a single generator, to reach the required accuracy. All these developments 
are by now under consideration. 
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